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We present a new analytical framework for describing the dynamics of a gravitational binary 
system with unequal masses moving with arbitrary relative velocity, taking into account the back- 
reaction from both compact objects in the form of tidal deformation, gravitational waves and self 
forces. Allowing all dynamical variables to interact with each other in a self-consistent manner 
this formalism ensures that all the dynamical quantities involved are conserved on the background 
spacetime and obey the gauge invariance under general coordinate transformations that preserve 
the background geometry. Because it is based on a generalized perturbation theory and the im- 
portant new emphasis is on the self-consistency of all the dynamical variables involved we call it a 
gravitational perturbation theory with self-consistent backreaction (GP-SCB). 

As an illustration of how this formalism is implemented we construct perturbatively a self- 
consistent set of equations of motion for an inspiraling gravitational binary, which does not require 
extra assumptions such as slow motion, weak-field or small mass ratio for its formulation. This case 
should encompass the inspiral and possibly the plunge and merger phases of binaries with otherwise 
general parameters (e.g., mass ratio and relative velocity) though more investigation is needed to 
substantiate it. 

In the second part, we discuss how the mass ratio can be treated as a perturbation parameter in the 
post-Newtonian effective field theory (PN-EFT) approach, thus extending the work of Goldberger 
and Rothstein for equal mass binaries to variable mass ratios 77 = m/M {M > m). As the mass 
ratio ri decreases from 1 to near 0, the smaller mass can be accelerated to higher speeds thereby 
requiring higher order post-Newtonian corrections to describe the system with sufficient accuracy 
for the data analysis requirements of current and future gravitational wave interferometer detectors. 
We provide rough estimates for the higher post-Newtonian orders needed to determine the number 
of gravitational wave cycles, with a specified precision, that fall into a detector's bandwidth. For 
example, we find that the number of cycles may need to be calculated through 4.5PN for ~ 5% 
precision (relative to 2.5PN corrections) for a compact binary with IOMq and lOOAf© constituents. 



I. INTRODUCTION 

Ground-based gravitational wave interferometer detectors, such as LIGO [1], are currently operating at their in- 
tended design sensitivity and are actively constraining the amplitudes of gravitational waves that may be passing 
within their detectable bandwidth. These observatories are expected to detect gravitational wave signals from the in- 
spiral, plunge, merger and ringdown phases of binaries with astrophysical masses, ranging from the mass of a neutron 
star to the mass of several hundred solar masses. Space-based detectors, such as LISA [2], which are currently in the 
planning and development stages are designed to detect the gravitational waves from inspirals of extreme mass ratio 
binaries and from binaries with nearly equal and large masses (e.g., with masses ~ IQ^Mq), among other sources. 

A. Existing analytical formalisms for gravitational binaries 

Different descriptions of gravitational binary systems ultimately depend on the parameters of the system. For 
example, the post-Newtonian (PN) approximation (see below) is useful for binaries with nearly equal masses (^ XQMq) 
during the inspiral phase, with extreme mass ratios where the small compact object moves in the weak-field region of 
the supermassive black hole, and with supermassive black holes orbiting slowly in a weak-field regime. For motion in 
a strong field, one often uses gravitational perturbation (GP) theory on a fixed but curved background (see below) 
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to describe the perturbations of a small compact object moving in the background field of a supermassive black hole. 
We review these formalisms in somewhat more detail below. 



1. The post-Newtonian expansion 

The post-Newtonian approximation is based on the assumption that two weakly gravitating objects orbit about each 
other at nonrelativistic speeds on a flat background. The strict weak field requirement can be lifted by constructing 
an annulus about each compact object where the PN-expanded metric is matched onto the near-field perturbed metric 
of the strongly gravitating compact object [3l, ^ . 

Iteratively solving for the metric perturbations and the positions of the masses yields approximate expressions in 
powers of the relative velocity for the phase of the emitted gravitational radiation, the energy and angular momentum 
they carry, the location of the innermost stable circular orbit, etc. To date, the equations of motion for the compact 
objects have been computed to order v"^ beyond the Newtonian ones, also denoted as 3.5PN, for non-spinning compact 
objects. For good reviews of the FN expansion see and [7]. 

The FN expansion is valid in a region, called the near zone, smaller than the wavelength of the emitted gravitational 
waves. Applying the FN approximation far from the source (in the far zone) can lead to logarithmic divergences since 
retardation effects from the finite propagation speed of gravitational waves cannot be neglected far from the sources. 
As a result, the FN metric in the near zone is matched to a metric containing the gravitational radiation emitted 
by the binary. This matching is performed in an overlapping region of the near and far zones. In this way, one 
can describe the generation and propagation of gravitational waves by slowly moving compact objects in the FN 
framework [S] H] . 

The metric in the far zone can be calculated from a post-Minkowski (FM) expansion where the Einstein equations 
are expanded in powers of Newton's constant Gn and solved iteratively for each order of the metric perturbations. In 
particular, there is no constraint on the velocities of the sources. This method is valid in those regions of spacetime 
that are weakly gravitating so the FM expansion is not applicable for compact binaries everywhere in the spacetime. 
However, as discussed above, the FM expansion can be applied far from such a system by matching onto the perturbed 
metric of the near zone computed in the FN approximation. 

The traditional approaches to the post-Newtonian expansion outlined above have been efficiently streamlined by 
Goldberger and Rothstein 8 who borrow time- honored techniques from quantum field theory. For example, Feynman 
diagrams are utilized to construct the FN expansion order by order and the divergences typically encountered in 
a theory of interacting fields and point particles are conveniently regularized using dimensional regularization [5]. 
Furthermore, a systematic formalism is introduced in [8 that allows for a point particle description of an extended body 
by utilizing an effective field theory paradigm and including all possible extra terms (so-called non-minimal couplings) 
into the point particle Lagrangian that are consistent with general coordinate invariance and reparamaterization 
invariance. 

Recent years have seen the development of the so-called effective one-body (EOB) approach [10] HI] that uses 
the 3FN expanded Hamiltonian to resum the conservative terms by mapping the two-body dynamics onto a system 
describing a test particle moving in an effective background spacetime. This background is a deformed Schwarzschild 
spacetime where the deformation parameter is the symmetric mass ratio, v = m M/{m + M)2. By following the 
geodesic motion of the test particle on this effective geometry one can describe the conservative dynamics of the 
binary even in the relativistic regime. The effects of radiation reaction are "grafted" onto the conservative test 
particle dynamics by using a flux balance argument where the power emitted in gravitational waves equals the loss 
in the particle's mechanical energy due to the reactive forces that come from emitting the gravitational radiation. 

While the EOB method has proved remarkably successful for providing approximate analytical waveforms of nearly 
equal mass non-spinning binaries it has several shortcomings from both a practical and theoretical point of view. On 
the practical side, the EOB approach does not yet describe well the dynamics of spinning binaries. Furthermore, for 
non-spinning binaries, one must patch the inspiral/plunge phases to the merger and the merger to the ringdown phase 
because the EOB approach does not smoothly describe these transitions. On the theoretical side, the EOB approach 
can be regarded as a phenomenological model with parameters that are obtained by fitting to numerical simulations 
of inspiral waveforms, for example. See [T21II31II1] and references therein for further discussions. 

2. Black hole perturbation theory for extreme mass ratios 

The extreme mass ratio inspiral (EMRI) scenario consists of two bodies with largely disparate masses. One body, 
a supermassive black hole (SMBH), has a mass so much larger than the other, a small compact object, that the 
dominant geometry of the spacetime is determined by the SMBH. Despite having a very much smaller mass than 
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the first, the small compact object nevertheless perturbs the background black hole spacetime and brings about 
the emission of gravitational radiation which causes the smaller mass to eventually spiral in toward the SMBH. Its 
backreaction on the small compact object generates a radiation reaction force which is called the 'self force'. It is a 
result of gravitational waves back-scattering off of the background spacetime and encountering the compact object at 
a later time and place in its orbit. This intrinsically nonlocal (dependent on the object's pass history) property is the 
main challenge in the calculation of the self force in this class of problems. See [15] for a good review of this approach. 

The perturbation theory for EMRIs has the advantage of treating the relativistic motion of the small compact object 
in a strongly curved spacetime. For detecting gravitational waves the space-based gravitational wave observatory LISA 
only requires knowing the self- force and the radiation through first order in the small mass ratio. For precise parameter 
estimation the self-force and the gravitational radiation will likely need to be calculated to second and third orders, 
respectively [IB] . 

B. New analytic methods for arbitrary mass ratios and relative velocity 

The approaches briefly reviewed above are largely useful for describing either nearly equal mass binaries (PN) or 
extreme mass ratio systems. While these constitute a large class of expected observable gravitational wave signals 
for both ground-based and space-based interferometer detectors there may also be signals coming from compact 
binaries having generally unequal or intermediate mass ratios. Even though the event rates of such systems are not 
accurately known, intermediate mass ratio inspirals (IMRIs) may occur in the dense centers of globular clusters |17j . 
Furthermore, during the final stages of galaxy mergers there may potentially be detectable gravitational wave signals 
from the inspiral of generally unequal SMBHs that originally formed the cores of the merging galaxies. Indeed, the 
potential discovery of an intermediate mass ratio binary [T5] (see also [T^ for an alternative interpretation of the 
data) suggests that IMRI sources may be present throughout the universe. 

Not much work has been done on IMRI systems, not because it is deemed less important or urgent but usually 
considered a very difficult task, as it falls in a hard- to- reach territory between the two familiar terrains, namely, the 
PN expansion and black hole perturbation theory. Simulating the dynamics of IMRI systems is not yet possible with 
numerical methods. As the mass ratio decreases, the computational cost and incurred numerical errors both increase 
dramatically making it difficult and very expensive to numerically evolve binaries with relatively different masses. For 
this reason even with small gains it would still be worthwhile to explore new analytic methods for inspiraling binaries 
having general mass ratios and relative velocities. Here we present two new analytical techniques to two different 
parameter regimes of inspiraling binaries with unequal masses, including IMRIs. 

1. Gravitational perturbation theory with self-consistent backreaction (GP-SCB) 

We have at least two motivations for inventing a new formalism. The first comes from a practical aim: to provide 
a common framework encompassing these two major schemes - the (equal mass) post-Newtonian scheme and the 
(extreme mass ratio) perturbation theory, but extending the limitations from both ends (of the mass ratio parameter) . 
We want this framework to be applicable to intermediate mass ratio processes without the slow motion, weak field 
restrictions. 

The second motivation comes from a more idealistic goal: to build a fully self-consistent formalism of perturbation 
theory that is able to account for the backreaction from all of the dynamical variables and on all of them. These 
variables (in a perturbative context here) are the motion of the compact objects (which is affected by radiation 
reaction and self- force from gravitational wave emission and by the tidal forces each object exerts on the other), the 
gravitational waves emitted by the binary, and the background metric (which together with the gravitational waves 
give the full metric). The full metric is dynamically determined by the negotiation between the two compact objects 
through mediation by the gravitational waves. 

To aid in such an effort we describe a generic framework for generating perturbation theories that is self-consistent 
in the sense that the participating dynamical quantities are conserved (in each perturbative order) on the background 
spacetime and have the appropriate behavior under general (i.e., not necessarily infinitesimal or small) coordinate 
transformations. This framework is based on previous work of Anderson |20j for studying the self-consistency of 
the gravitational geon solution in [21]. As a useful check of this new formalism we demonstrate how to obtain the 
well-established post-Newtonian expansion and the extreme mass ratio perturbation theory as its subcases. 

To apply this new formalism it helps if one can identify some physical parameter in certain stages of the binary's 
motion with a clear discrepancy in the mass, velocity or frequency most suitable for a perturbative expansion |41j . As 
an example, we consider those cases where the dynamics of the binary includes the existence of an inspiral phase. An 
inspiraling binary evolves with an orbital period that is shorter than the secular time scale associated with the radiation 



4 



reaction and dissipative self-forces. We demonstrate that averaging over the fast orbital period leads to a perturbation 
theory where the mutual backreaction from and on all of the dynamical variables involved are incorporated. In this 
case the companion-induced tidal moments are included into the description of the compact object (s) [5]. Because 
this perturbation theory does not depend on the binary's mass ratio or relative velocity, at least explicitly for its 
formal construction, we have reasons to believe that our treatment in this example may be applicable to inspiraling 
binaries with generic masses and velocities, including intermediate mass ratio binaries and comparable mass binaries, 
and possibly to the plunge and merger phases because of the inclusion of backreaction effects. Additional assumptions 
and approximations may be useful for practical calculations but they are not necessary for the internal consistency 
of this new theory. 

The self-consistency requirement is a crucial structural feature of this new approach, and the allowance for dynamical 
backreaction is an attractive functional feature because the background can respond to the effective stresses and 
energies arising from the motion of the compact objects and the emitted gravitational waves. This is to be contrasted 
with other prevailing approaches, foremost the PN and EMRI perturbative schemes, that choose a fixed background 
which never deviates from its originally specified form. While this is clearly convenient for calculations of many 
scenarios, especially if the fixed background possesses some isometrics, it is not general enough for a wider class of 
binary processes and not very realistic on a broader scope, because of the nonlinear interactions between the massive 
objects and the nonlocal (history dependent) influences of the gravitational waves. 



2. Post- Newtonian effective field theory (PN-EFT) for unequal masses 



The PN expansion is often regarded as being applicable to nearly equal mass binaries or to the weak-field evolution 
of an extreme mass ratio binary. However, many interesting gravitational wave signals for LIGO and LISA come from 
evolution in a strong field regime, including the relativistic motion of a compact object near a supermassive black 
hole and the plunge and merger phases of an astrophysical equal mass binary. 

To develop equations of motion that are based on a PN approximation but can also be applied to unequal mass 
binaries in stronger fields requires introducing the mass ratio r] = m/M as an additional expansion parameter to 
the velocity v and developing the current PN expressions to higher orders. Within the context of the effective field 
theory approach originally developed by Goldberger and Rothstcin "5", who called this approach 'non-relativistic 
general relativity' (NRGR) [42J, we can include rj into the power counting scheme for determining the scaling of the 
interaction terms that contribute to the PN expansion at a given order in v and r/. In fact, in the context of PN-EFT, 
this can be done systematically and efficiently to any order. 

For a given mass ratio and total mass it is instructive to have some idea of the PN order that a quantity, such 
as the number of cycles of a gravitational waveform that fall into LIGO's bandwidth, should be expanded through. 
With such estimates one can begin to more accurately compute PN expanded equations of motion, waveforms, etc. 
for binaries that have mass ratios of 10~^ or 10~^, for example. For such systems, the smaller compact object can be 
accelerated to higher speeds and so a knowledge of higher order PN terms becomes important for accurately describing 
the binary. 



C. Organization 

This paper is organized as follows. In Section [IT] we present this new perturbation theory for gravitational binary 
systems with self-consistent backreaction from and on all dynamical variables involved within a rather general frame- 
work, viz., the compact object motion, its tidal deformations, gravitational waves, and the background spacetime 
geometry. We demonstrate how the post-Newtonian expansion and the perturbation theory for describing EMRIs 
arise from within this formalism. As an example, we apply this new formalism to dynamics which has an inspiral 



phase for any mass ratio and relative velocity. In Section III we generalize the PN-EFT approach to include the mass 



ratio as an expansion parameter. We then provide crude estimates that indicate the post-Newtonian orders needed, 
with a given precision, to calculate the number of cycles observed in LIGO's bandwidth for a binary with general 
values for its masses. We end with a discussion of directions for future development. Further details for each of these 
approaches are contained in two forthcoming papers [221 123] . 



II. GRAVITATIONAL PERTURBATION THEORY WITH SELF-CONSISTENT BACKREACTION 



Our goal is to establish a formalism which describes the evolution of three dynamical variables in a self-consistent 
manner: the spacetime and its perturbations (the gravitational waves), and the motion of the compact objects. We 
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take a first principles approach that is sufficiently flexible and general to accommodate several approximation schemes 
and perturbation theories. This general theory should encompass all the existing approaches as delineated in the last 
section yet be able to treat parameter ranges outside of their validity. There are a few basic criteria we set for any 
such formalism to obey, foremost a gauge-invariant effective stress-energy tensor for gravitational waves. For these 
purposes we adopt the framework of Anderson [20j who provided a rigorous foundation for the Brill-Hartle-Isaacson 
averaging procedure [Ml 123 ESI and apphed it to a careful study of gravitational geons [21] . 

The two-body problem may be described at the level of the equations of motion in several ways, depending on the 
physical setup, the processes, and what specific parameter values are of interest. There are three classes of compact 
binaries (i.e., radiating binary systems composed of compact objects) typically regarded in the literature: 1) black 
hole/black hole (BH/BH), 2) black hole/neutron star (BH/NS), and 3) neutron star/neutron star (NS/NS). 

In a first principles description, the spacetime geometry of a BH/BH binary is described by the vacuum Einstein 
equations 

GM = (2.1) 

where g^i/ is the full metric of the spacetime. When a binary is composed of one or two neutron stars one can describe 
the system in terms of a matter stress-energy via the non-vacuum Einstein equations, 

G,,.(.9) - ^^G^Tl^fig- p,p, . . .) (2.2) 

where the summation is included depending on whether one is considering BH/NS or NS/NS binaries. Here g^^i, is 
the full spacetime metric, which is determined by the dynamics of the binary, and T^'^ is the stress-energy tensor 
appropriate for a neutron star with density p, pressure p, etc. At present, the dependence of this stress tensor on the 
fluid variables is not well understood because the NS equation of state is not known. 



Solving (2.1 1 or (2.2) can only be done using numerical techniques. Even then, most simulations are carried out for 
nearly equal mass binary black holes (i.e., mass ratios > 1/4) and numerical methods for BH/NS and NS/NS binaries 
are still in their infancy. 



Analytically solving (2.11 or (2.2) can be accomplished by invoking approximations and restricting solutions to a 
particular dynamical regime (e.g., inspiral or ringdown phases). In a BH/BH binary having one mass considerably 
smaller than the other, as in the extreme mass ratio scenario, one can model the smaller BH using a point particle 
description. In this case, the Einstein equation ( |2.1[ ) is approximately 

G^,{g)^8nGNTPP{g;z) (2.3) 

where the coordinates of the particle worldline are denoted by z^. Here is (close to but not quite exactly) the full 



metric of the binary spacetime (as obtained from 2.1). For very small mass ratios one could provide an equivalent 
description to a high degree of accuracy by replacing the binary system by a background spacetime generated by the 
larger BH otherwise in isolation plus a perturbation generated by the presence and motion of the smaller BH. As 
such, one may decompose the full metric into a flxed background go/ji^ (e-g-, describing a Kerr black hole) and its 



perturbations h^j^^. By expanding (2.3) in powers of /i^i^ and the mass ratio one can perturbatively obtain analytic 



solutions for a binary BH in the extreme mass ratio limit. 

From the considerations above, the equations of motion for a general compact binary system can be generically 
written as 

G^A9,(^) = (2.4) 

where 

Q^iv{g,a) = G^u{g)-%„{g,a) (2.5) 
and 7^1/ is proportional to the stress-energy for the compact object(s), if the assumptions and input allows for the 



inclusion of material stress-energy. In (2.5), a denotes the collection of variables that describe the dynamics of the 
compact object(s), such as the density, pressure, worldline coordinates, etc. The number of contributing stress tensors 
to T^n and the forms of these contributions depend crucially on how one describes the constituents of the binary (e.g., 
as a point particle, as an extended fluid body), the regime of dynamical evolution being considered (e.g., inspiral, 
ringdown), etc. Below we closely follow the presentation given in [ 20] . 

A. Generalized gauge transformations, gauge invariance, and self-consistency 

In an analytic approach the metric of the full spacetime g^i, is separated into a background part g^^u and a 
perturbation part /i^j/ so that g^j^^ = (?o^ii^ + We stress that this decomposition is completely arbitrary. Indeed, 
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in standard perturbation expansions one fixes the background metric according to the problem at hand (e.g. Kerr 
background for EMRIs). The unperturbed Einstein equation (2.4 1 can generally be written according to the above 
metric decomposition in the form |43j 



or equivalently as 



G,,uigo, a) + AC/^^(/i, a) 



G{go) + AG{h) = T{go, a) + AT{h, a) 



(2.6) 



(2.7) 



Because we want the metric decomposition to be arbitrary the quantities AG^jy and AT^^ contain all of the dependence 
on the gravitational perturbations /i^^ and are not necessarily small with respect to the background Einstein or stress 
tensors, G^i,{gQ) and T^,^{go,a), respectively. 

For the formalism to be self-consistent we require that the perturbed quantity 



Ag{h, a) = AG{h) - AT{h, a) 



(2.8) 



is both conserved and invariant under coordinate transformations that preserve the structure of the background 
geometry. The former is straightforward to demonstrate. 

Given a solution {gof^^^a) to (2.7 1 we know that G(.go) and T{gQ^a) are separately conserved with respect to the 
background geometry upon using the Bianchi identities for the Einstein tensor and the conservation equation for the 
stress tensor, which gives the equations of motion for the compact object variables on the background geometry. It 



therefore follows from (2.7l that (2.8l is also conserved with respect to the background geometry, V^A5^^(/i, a) = 



where is the covariant derivative compatible with 50^1^- 

Since the perturbed quantity AQ(h, a) is not necessarily small with respect to Q{go) it follows that the coordinate 
transformations, which change the perturbed metric but not the background, are not necessarily infinitesimal as is 
usually considered for discussions of gauge invariance. As such these coordinate transformations are called generalized 
gauge transformations [20]. 

A coordinate transformation can always be written in the following form 



(2.9) 



where is a function of and is not necessarily infinitesimal nor small. Accordingly, the full metric changes via 
the usual tensor transformation rule so that the metric perturbations in the two coordinates are related through 



+ [ffOa/5(x)+/la/5(5)]e" e''. 



(2.10) 



where the derivatives of are with respect to x^ . In the limit of infinitesimal we recover the usual form for the 
infinitesimal coordinate transformation of the metric perturbations 



hf_,y{x) = h^y{x) - gaf_,„M{x)C - 9op.a{x)^"^ - goau{x)^° 



(2.11) 



The perturbed quantity AQ{h,a) is gauge invariant under a generalized gauge transformation (2.9 1 if the following 
relation holds 



Agih,a) = Ag{h,a) 



(2.12) 



where /i^i^ is found by inverting (2.101. The proof can be found in '20]. 

So far we have described a general framework that is self-consistent in the sense that AQ is both conserved on the 
background geometry and is gauge-invariant under coordinate transformations of the form (2.9 1. However, we cannot 
solve for the background or the perturbations from the Einstein equations 



^/(5o,a) = -AQ{h,a) 



(2.13) 



alone without knowing how the full spacetime geometry is decomposed into a background and its perturbation. There 
are at least three ways to do this, in principle. The first is to specify the background and solve for the perturbation 
or specify the perturbation and solve for the background. 

The second is to define a conserved and gauge invariant effective stress energy tensor (or in some approximate 
sense) for gravitational waves T^™ on the background so that 



Ag{h,a) 



-8TrGNTS'"{h,a) 
8TrGNTS'^{h,a) 



(2.14) 
(2.15) 
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For example, one could define 



\ttG 



N 



-{Ag{h,a)) 



which is clearly gauge invariant under coordinate transformations of the form (2.9 1 
a time, space, or space-time average. We note that this stress tensor is related to t 
in |24j with time averaging and by Isaacson j25l I26j with space-time averaging. 



(2.16) 

Here the angled brackets denote 
le choice made by Brill and Hartle 



The third is to introduce a gauge-invariant equation that determines the metric perturbation thereby fixing the 
spacetime decomposition. For gravitational waves this is naturally provided by a wave equation 



Higo,h,a)^0 



(2.17) 



We will discuss this approach in more detail below with a couple of examples. 

In most situations encountered in the two-body problem the equations cannot be solved exactly and some ap- 
proximations must be invoked, usually in the form of perturbation expansions. Let us assume that there exists an 
expansion of AQ so that 



Ag==Aig + A2g 



(2.18) 



Here the subscript on a A denotes the order of the expansion. We remark that we have not actually specified any 
particular assumptions regarding the expansion such as the strength of the gravitational waves or its derivatives nor 
have we identified a small parameter to do the expansion. We merely assume that there exists a valid expansion of 
the form (2.181. As such, it is not necessarily true that A2G is strictly quadratic in h, for example, as is the case for 



the gravitational geon |20l [25 • 

If we truncate the expansion in (2.181 to order n then the discussions abo ve imp ly that AiQ + - ■ ■ + AnG is conserved 
with respect to the background geometry to n"^ order. Gauge invariance of (2.18) can also be shown provided that the 



generalized gauge transformations are such that h and h are of the same order of magnitude to facilitate comparison 
[5ni[lSl[53], which we will assume from here on. It is important to realize that the quantity AiQ + ■ ■ ■ + AnG is gauge 
invariant to n*'' order. Therefore, AiG is gauge invariant to first order but A2G is generally not gauge invariant at 
any order. Instead it is the sum AiG + A2G that is gauge invariant to second order. 



B. Existing theories as special cases 

In this section we will show that the general perturbative framework developed above contains the familiar existing 
theories including 1) the (Regge- Wheeler) perturbation theory for EMRIs and 2) the post-Newtonian expansion for 
equal mass binaries. In the following section we will generate a new expansion that may be largely insensitive to the 
mass ratio and the typical velocities of the binary. 



1 . Perturbation theory for EMRIs 

To build a perturbation theory for EMRIs, use the smallness of the binary's mass ratio e = m/M <C 1 so that 
hfj,i, — 0{e) and TPP{go, z) — 0{e) where T^^ is the stress tensor of a point particle, which models the small compact 
object. Notice that we are not saying anything yet regarding the nature of the background spacetime, namely, whether 
it is fixed or dynamical, etc. Expand the Einstein equations after decomposing the full metric as g ^ go + h to find 

Gigo) + AG{h) = AT{h, z) = 87rGArTfP(5o, z) + 87rGArATPP(/i, z) (2.19) 

Notice that the source of curvature should be regarded as a perturbed quantity since the point particle stress tensor is 
proportional to the small mass m of the compact object. Expanding the perturbed quantities AG and AT in powers 
of e gives to second order 

AG(/i) = AiG(/i)+A2G(/i)-h--- (2.20) 
Ar(/i, z) = AiT{z) + A2T{h, z) + • • • (2.21) 

where Ai^2 denotes a quantity of 0(e^'^). Because the particle stress tensor is already proportional to m then AiT(z) 
is independent of h. 
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We can fix the metric decomposition by specifying a gauge invariant but otherwise arbitrary equation for the 
perturbations, H(go, h, z) = 0. A natural choice for iJ is a wave equation such that 



H{go, h, z) = AG{h) - AT{h, z) = 



From (2.22) and (2.191 it follows that the leading order contribution to the Einstein equations is 

G{go) = 



(2.22) 



(2.23) 



Note that according to ( |2.23 1 the background geometry is not affected by the gravitational waves. Indeed, the 
background metric describes a vacuum spacetime that one can fix or specify ab initio. To solve the wave equation 
order by order in e we write 



so that (2.221 becomes 



h=hi + h2 + 0{e^) 



AiG(/ii) = Air(z) 
AiG(/i2) = -A2G(/ii) + A2T{hi,z) 



(2.24) 



(2.25) 
(2.26) 



through second order in e. The first line simply describes the wave equation for the first order perturbation hi, which 
is sourced by the stress tensor for the point particle with worldline coordinates z^. 

The particle's equations of motion follow from the geodesic equation 0^(17, z) — due to the conservation of AT on 
the full spacetime g. Using the metric decomposition g = go + h and the expansion of h in (2.24) we see that through 
first order 



a{go,z) 

which can be written in a more recognizable form as 



-Aia{hi,z), 



u"(t)V„u''(t) = - Vo.hi'^Jz) - -V^hi^piz) M"(r)u'^(T) 



(2.27) 



(2.28) 



where r is the proper time for the worldline 2^ and — dz'^/dr. Indeed, this is the unregularized equation describing 
the first order self- force on the small compact object, which was first derived (and regularized) in [27l [28]. 



2. Post-Minkowski and post- Newtonian expansions 

A post-Minkowski expansion is built by including a stress tensor for each of the two compact objects and choosing 
Newton's constant Gfq as the expansion parameter [ll]. We find the following self-consistent perturbed field equations 
in a manner similar to the previous section 

G(.go) = (2.29) 
AiG{hi) = AiT{z,Z) (2.30) 
AiG(/i2) = -A^Gihi) + A2r{hi,z, Z) (2.31) 

where z^ and denote the worldline coordinates of the masses m and M , respectively, and 

Air(z, Z) = 87rG^T(''^)(5o, z) + %^GnTII,^{9^, Z) (2.32) 
A2T(;ii,z,Z) =87rGwAiT^^)(/ii,z) + 8^G^AiT^,)(;ii,Z) (2.33) 

Again, the leading order equations indicate that the background is vacuous and the considerations of the specific 
problem imply g^^^ = 77^^. 

In the near-zone region of a comparable mass binary, where the gravitational perturbations propagate nearly 
instantaneously, the compact objects move slowly compared to the speed of light. Consequently, we use the relative 
velocity v as an expansion parameter. Together with the post-Minkowski expansion (where we assume that Gpf — 
0{v^) from the virial theorem) we can construct a valid perturbation theory of the Einstein equations and particle 
equations of motion. 
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Following the similar steps as in the previous section we find at leading order that G{go) = 0, which indicates that 
the background can be fixed to Minkowski spacetime as before and that the gauge invariant wave equation is 

AG{h,z) = AT{h,z,Z) (2.34) 

At first order in v we find the perturbed wave equation is AiG(/ii) = where the right hand side is zero because the 
leading order contribution from each point particle stress tensor is 0{Gn) = 0{v^), which is a second order quantity. 
Therefore, we can consistently choose hi = 0. At second order in v the perturbed wave equation is 

AiG(/i2) - A2Tih, z, Z) - 8^GArT(„,)(77, z) + 8nGNT(^M){v, Z) (2.35) 

where we have used that hi = 0. One can easily show that when a gauge is chosen for the metric perturbations (e.g., 



the harmonic gauge on a fiat background) that (2.351 is just the Poisson equation for /iqo sourced by the two point 
particles on worldlines with coordinates 2'^ and Z^^. 

The particle equations of motion follow from the geodesic equation in the full spacetime so that through 0{v^) we 
find 

a{m){v, z) = -Aia(„)(/i2, z) (2.36) 
a{M){v, Z) = -Aia(M)(^2, Z) (2.37) 

which are the equations for two particles moving in their respective Newtonian gravitational potentials. Continuing 
to higher orders merely reproduces the equations of motion for the field and the particles at successively higher 
post-Newtonian orders. 

C. Two-time separation and adiabatic expansion in the GP-SCB formalism applied to the inspiral phase 

In this section we illustrate how the general perturbation formalism with self-consistent backreaction can be applied 
to a class of problems with the introduction of a two-time separation scheme. Throughout this section we consider the 
motion of a compact object with mass m (a neutron star or black hole) in the spacetime of a larger black hole with 
mass M > m. While it is convenient to use a point particle description for the smaller compact object it is difficult to 
justify such an approximation for an extended body. Even if it is acceptable for a certain period of time one cannot 
assume this is valid throughout the entire course of its evolution for binaries with arbitrary mass ratios because the 
larger mass will induce tidal effects that manifest as higher order multipole moments for the distorted shape of the 
smaller compact object. Such effects are not negligible during the plunge and merger phases, especially as m —^ M. 

To include the tidal effects in such binaries one needs to take into account the backreaction effects. For comparable 
mass binaries the compact objects are more severely disrupted during the inspiral, plunge and/or merger stages in 
which case a point particle (pp) description for the masses is more likely to become invalid. This is where an effective 
point particle (epp) description beyond the pp approximation will prove useful in the GP-SCB scheme. 

The non-minimal couplings that appear in the effective point particle description of an extended body [RJ are related 
to terms in a multipole expansion that scale as powers of rm/A where is the size of the compact object with mass 
771 and A is the typical wavelength of an external gravitational perturbation at some scale. For such a compact object 
plunging toward merger with another mass M, 

"jlL ^'J^y<^'j:^ (2.38) 
A r r 

for arbitrary velocities (The second relation comes from the fact that for both EMRI and equal mass binaries, 
V « 0.5 — 1.0 near merger). Near the merger phase the orbital separation is of order r ^ r„i + rj\/ and it follows that 

Thus if the mass ratio is <^ 1 then one can justify using a point particle description for the smaller mass so long as 
a sufficient number of non-minimal terms (possibly many) are included to account for the companion-induced tidal 
moments. 

The action describing the effective point particle (epp) dynamics for the small compact object is given by [5] 
Sepp[z,g]= — m J dr + cji J dr R{z) + cy J dTRap{z)z"z^ 

+ CE j drSM^'^^i-A + CB I dTB^f,{z)B"''\z) + --- (2.40) 
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where the quantities £ap and Bap are the electric and magnetic parts, respectively, of the Weyl curvature tensor. The 
constants {ci\ depend on the internal structure of the compact object (e.g., parameters of the equation of state for a 
neutron star) and can be calculated by matching onto a full description of the "microscopic" theory for the compact 
object. See [5] for further discussion. From this effective point particle action one can obtain the stress tensor for the 
smaller compact object 



from which we obtain the Einstein equation for the binary 



(2.41) 



(2.42) 



Because of the non-minimal couplings in (2.40) the motion of the compact object to is not given by the geodesic 



equation. In fact, the acceleration of the smaller mass to moving in the full spacetime is given by the conservation of 
j^epp ^Yie full spacetime. 



(to - crR{z) + • • • )a^(r) = cii{g^'' + i^'i-^) V,i?(z) + 



(2.43) 



and is non-zero precisely because the compact object is an extended body although it is being modeled with a point 
particle description. Notice the appearance of an effective mass resulting from these non- minimal interactions |45j . 
We can write the particle equations of motion (2.43) in a general form 



fn^^^g, z)ai,{T) = i^^(g, z; cr, cv , ■ ■ ■) 



(2.44) 



where the effective mass is a tensor that is generally time and space dependent, and F^^ accounts for the forces on the 
particle arising from the finite size effects of the tidally distorted compact object m. When such companion-induced 
tidal moments can be ignored (2.44) reduces to the geodesic equation on the full spacetime with metric g. 

Binaries exhibiting an inspiraling phase admit the characteristic feature wherein the orbital period Torb is often 
much shorter than the timescale for radiation reaction and dissipative self- force effects Trr. Therefore, let us introduce 
an averaging procedure that effectively averages the fast motion of the binary. 

Split the spacetime into a background and its perturbation so that g = go + h, and write the worldline coordinates 
of the effective point particle as z — zq + where zq denotes the coordinates of the "background" worldline and its 
perturbations. Then, introducing e as the expansion parameter, defined as the ratio Tort/TuR, which is small during 
at least the inspiral phase, we can expand (2.42) and (2.44) in a series giving 



G(5o) + AG{h) = SnT-PPigo, zo) + SttAT^p^ (/i, C) 
m{go, zo)a{go, zq) + m{gQ, zo)Aa{h, C) + Am{h, C)a{go, zq) + Am{h, C)Aa{h, () = F{gQ, zq) + AF{h, () 

where we have dropped the spacetime indices in the particle equations of motion for notational simplicity. 
Expand the perturbed quantities AG and AT^pp to second order 

AG ^ AiG + A2G + 0{e^) 
AT^PP = AiT^PP + A2T''PP + 0(6^). 



(2.45) 
(2.46) 



(2.47) 
(2.48) 



We fix the metric decomposition by imposing a condition on the metric perturbations that uses an averaging procedure 

(AiG(/i) - 87rGArAir"PP(/i, 0) = 0. (2.49) 

We have in mind that the average is over an orbit of the evolution. In particular, the average should be defined so 
that the secular decrease in the orbital period is incorporated in its definition. However, the precise definition of the 
average used here is not necessary for formulating the perturbation theory. 

Notice that the condition in (2.49) is not exactly gauge invariant, which is simply a manifestation of the fact 



that fixing any metric decomposition is an intrinsically gauge-dependent procedure. As such, (2.49) corresponds to 



choosing a particular gauge and so will not necessarily remain invariant under a generalized gauge transformation of 
the form (2.9 1. Nevertheless, the quantity 



(AiG(/i) + A2G{h) - 8TrGNAiT'PP{h, () - 87rGNA2T^PP{h, ()) 



(2.50) 



will be gauge invariant through second order, as mentioned in earlier discussions. The condition (2.49 1 is a reasonable 
choice since the fast and slow time scales of the orbital motion are imprinted into the gravitational variables. In a 
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sense, we are assuming that the slow variation of the gravitational field (~ Tra) is contained in the background go 
while the relatively fast variations (~ Torb) are contained in the metric perturbations h. 
Further expanding the metric perturbation as 



h=hi + h2 + O(e^) 



(2.51) 



gives for (2.49) to second order 



AiG(/i2) - A2G(/ii) (2.52) 
AirPP{h2,C) - A2T'=fP(/ii,C) (2.53) 
where ~ is to be read "is of the same order as." Using ( 2.51[ ) we can expand (2.491 in e to find, order by order, 

(AiG(/^i) - 8TTGNAiT^PP{huO) = (2.54) 
(AiG(/^2) - 8TrGNAiT^PP{h2,0) = (2.55) 



the relation 



through O(e^). Notice that we use (2.491 to define the gravitational equations of motion in a particular coordinate 
system and use (2.51) to iteratively solve these in the same coordinates. 



Just as (2.49) defines the metric decomposition go + h we now define the worldline decomposition + C through 



(Aim(/i, C)a(ffo, zo) + m{9o, zo)Aia{h, C) - AiF{h, ()) = 



(2.56) 



so that the averaged first order perturbed acceleration, effective mass, and finite-size effect force vanishes. For pure 
point particle motion, (2.56) reduces to ( Aiajh , C)) = 0. 



Writing C — zi + Z2 + 0{e^), averaging (2.45) and (2.46), applying the constraints (2.54) and (2.55) to these along 



with the definition of the averaging procedure in (2.56), and subtracting from the unaveraged original expressions 
(2.45) and (2.46) yields the following field equations for the gravitational variables 



G{go) = T'PP{go, zo) - {A2G{h^)) + 87rGjv(A2r^PP(/ii) + A2T-pp{z^)) (2.57) 
AiG(/ii) = %TTGNAiT''PP{hi) + 8ttGnAiT''pp{zi) (2.58) 
AiG(/i2) = - A2G(/ii) + (A2G(/ii)) + B^GjvAiT^Pf (/i2) + 8ttGnAiT''pp{z2) + 87: GnAiT^pp {hi, zi) 



+ 8TrGN [A2T''PP{hi) + A2T'PP{zi) - (A2T^^'^'(/ii) + A2T^f^'(zi))] 
and for the worldline coordinates of the effective point particle 

m{go, zo)a{go, zo) = F{go, zo) - {A2m{hi) + A2m{zi))a{go, zq) - m{go, zo){A2a{hi) 
- ((Aim(/ii) + Aim(zi)) {Aia{hi) + Aia(zi))) + {A2F{hi) 4 
m{go, zo){Aia{hi) + Aia(zi)) = - (Aim(/ii) + Aim{zi))a{go, zo) + AiF{hi) + AiF(zi) 
m{go,zo)Aia{z2) = ■■■ 



(2.59) 

+ A2a(zi)) 

A2F(zi)) (2.60) 
(2.61) 
(2.62) 



(We have not written the second order equations for the particle coordinates Z2 due to their length.) 

The equations of motion ( 2.57 1-( 2.62) describe the inspiral phase of a binary. Since we did not build into the 



expansion any explicit reference to the mass ratio or the relative velocity these equations may be applicable for 
describing an inspiraling binary with arbitrary values for these quantities. Our hope is that these equations may also 
be useful for describing the plunge and possibly merger phases since we can incorporate the companion-induced tidal 
effects on the smaller compact object during the course of its evolution in the strong field region of its larger mass 
companion. While such a description using the effective point particle worldline of [5] will surely break down if the 
mass m undergoes tidal disruption or the horizons of two black holes just merge it may nevertheless be useful for a 
description as the merger phase is approached. However, further investigation is needed to understand these points 
and the domain of validity of this system of equations. 

The physics contained in these equations is very rich. The worldline coordinates, gravitational wave perturbations 
and background geometry mutually back-react on each other. This can be seen in (2.57), for example, where the 



background geometry is sourced by the leading order stress tensor for the effective worldline zq, by the averaged 
second order gravitational waves (analogous to an averaged Landau-Lifshitz pseudo-tensor) ^ (A2G(/ii)), and by the 
averaged higher order corrections of the particle stress tensor. As such, the background geometry cannot and ought 
not be specified but must be determined self-consistently with the dynamics of the other gravitational and particle 
variables. 

Notice also that the equations for zq are sourced by averaged forces involving the first order perturbed coordinates 
and gravitational waves. In particular, the worldline coordinates zq do not follow a geodesic of the background, even 
when the non-minimal couplings from (2.40) are removed, but are determined by consistently solving (2.60 1-(2.62 1. 
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III. POST-NEWTONIAN EXPANSION FOR ARBITRARY MASSES 



The effective field theory paradigm can be appHed to any gravitational system where the typical scales of the 
gravitational perturbations are much larger than their sources or scatterers. These conditions can be satisfied for 
compact binaries with an arbitrary mass ratio during some portion of its evolution (e.g., inspiral and ringdown 
phases). For many binary systems, such as a neutron star and a IOMq black hole, the mass ratio parameter rj is small 
suggesting that it can be used as an additional expansion parameter in the post-Newtonian approximation. In this 
section we discuss the role that the mass ratio 

TTl 

V=j^ , M>m (3.1) 

has in parameterizing the relevant interactions at the desired order in the effective field theory approach to the 
post-Newtonian expansion, PN-EFT. 

One purpose for doing this is to provide a better understanding for the description of binary systems when the mass 
ratio is tuned to smaller values. The analytical techniques for describing (nearly) equal mass binaries are considerably 
different than those used for binaries with extreme mass ratios. Whereas the former is based on assumptions of small 
velocities and weak gravitational fields the latter is grounded in black hole perturbation theory and is applicable 
to relativistic binaries with strong gravitational effects. Our approach, described below, can be regarded as a way 
to optimize the post-Newtonian expansion for binary systems with unequal masses by organizing the perturbative 
corrections into groups of terms that have numerically similar magnitudes and relevance at a given post-Newtonian 
order. 

Incorporating 77 into the approximation scheme allows for a more quantitative estimation of the magnitude of the 
interactions that contribute at a given PN order since the smallness of the mass ratio can suppress some interactions 
relative to others. For extreme mass ratio binaries (77 <C 1) the interactions accounting for the backreaction on SMBH 
are strongly suppressed relative to those describing the leading order motion of the small compact object, for example. 



A. Arbitrary mass ratios in PN-EFT 



Introducing 77 into the post-Newtonian expansion requires modifying the PN-EFT power counting rules originally 
developed for the (nearly) equal mass scenario in [|]. Let z{t) and Z(i) denote the trajectories of the masses m and 
M, respectively. In terms of these, the binary's center of mass R is described by 

R(t)^-^z{t) + -^Z{t). (3.2) 
1 + 7/ 1 + ?7 

In the center of mass frame the velocities of the two masses are related through 

V{t) = -7^v{t) (3.3) 

where v = z and V = Z are the 3- velocities of m and M , respectively. From the virial theorem of Newtonian gravity 

2 , 2GNmM 
m^^ + MV^ = , 3.4 
|z - Z| 

we find the following relations between the velocities of the masses and their respective leading order potentials, 

2 GjvM M 

V 2—, 3.5 

r 

2 Gf^m ^ M 

V =77 77 . (3.6) 

Here r denotes the magnitude of the orbital separation, r = z — Z, and we use the same conventions as [8]. In 
particular, we use units where c = h = 1 and m~i^ = 32^0^. 

In PN-EFT there are two kinds of metric perturbations: potential gravitons Hf^^, and radiation gravitons |46j . 
The orbital separation of the masses is much smaller than the light-crossing time so that each compact object interacts, 
to leading order, with the instantaneous potential iJ^jy generated by the other mass. At the orbital scale, the emission 
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of gravitational waves is described by an external, slowly varying, long wavelength perturbation. Therefore, the full 
metric is given by 

. Hfj^i, hfxu /„ _N 

g^iy = tlt^y H — + (3-7) 

m.pi rripi 

where the perturbations are canonically normalized by the Planck mass (equivalently, by Gjy ). 

The power counting rules are given by the virial theorem, (3.51 and ( |3.6| , and the ratios of the masses to the Planck 
mass, m/nipi = y/rjvL and M/nipi = yJvLji]. Here L ^ raw is the center of mass angular momentum of m and the 
scaling laws of the potential and radiation gravitons are the same as in 

The total action of the system is given by 

S'totiz, Z, }i\ = S\r] + H/rUpi + h/nipi] + S'epp[z, H, h] + Sepp[Z, iJ, h] (3.8) 

where the first term is the gauge-fixed Einstein- Hilbert action and the last two describe the masses (treated as effective 
point particles [8 ) coupled to potential and radiation gravitons. There are many non-minimal couplings in the effective 
point particle Lagrangians that describe the effects from companion induced tidal moments but these arc higher order 
effects that we will ignore for most of our discussion. 

The larger mass M couples locally to the potential graviton via factors of 

Hf,, {t, Z{t))^ H^, {t, -77z(i)) (3.9) 

in the point particle Lagrangian. The mass ratio parameter appears in the argument of H^^ and implies that all 
?i-point functions of potential gravitons will give rise to contributions proportional to products of 

(i/,.(M)i^.M^,Z))oc^^^^^, (3.10) 

which do not scale homogeneously with rj. However, by expanding the potential graviton in multipoles about the 
center of mass, which we take to be at the origin, 

H^,{t,Z{t)) = H^,{t,Q) + Z\t)d,H^,{t,Q) + ]^Z\t)Z^{t)d,d,H^,{t,Q) + ■■■ (3.11) 

we see that each term above now scales as a definite power of the mass ratio since {Z^ diY H ^ "rfH^i, for some 
non-negative integer p. As in j8] the radiation gravitons are also expanded in multipoles 

V = 0) + z'd,h^,{t, 0) + ^z'z^d,djh^,{t, 0) + • • • , (3.12) 

Z(t)) = V(t,0) + Z'd~h^,{t,0) + ^Z'Z^d,dfh^,{t,0) + ■■■ . (3.13) 

The power counting rules imply that z^dih vh whereas Z^dih rjvh. 

As an example of the new power counting rules, it is straightforward to show that companion-induced tidal effects 
(i.e. the effacement of internal structure) for the smaller mass ra first appear at order r/'^v^'^ and for the larger mass 
M at order rj^v^'^. In the EMRI limit, where v ^ 1, these agree with earlier results [50] . 



B. Estimates of post-Newtonian orders for the number of gravitational wave cycles 

Current and future generation gravitational wave interferometers will be capable of extraordinarily precise mea- 
surements partly because the interesting signals, which are buried in a sea of detector noise, can be extracted using a 
technique called matched filtering . Matched filtering is most efficient when the sought after signal is known 

or modeled with a sufficient accuracy. For gravitational wave interferometers the signal is the waveform 

h{t) = F+h+it) + F^hy,{t) (3.14) 

given in terms of the -I- and x polarizations of the gravitational wave and the pattern functions i^+,x, which depend 
on the geometry of the interferometer and the direction to the signal source. For inspiraling sources, these detectors 
will accumulate many oscillations of the gravitational wave so that the phase of h(t) must be known to a high degree 
of precision. 
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(a) 6 = 0.1 



{h)5: 



•1 

: 0.05 



(C)5: 



: 0.01 



FIG. 1: Contour plots showing the post-Newtonian orders q/2 needed to estimate iZiyc with a precision (relative to the 2.5PN 



contribution) of (a) 10%, (b) 5%, and|(c)| 1% as a function of the total mass (in Mq) and the mass ratio rj. The white areas 



indicate those binaries with masses that have ISCO frequencies lower than lOHz, which we regard here as not detectable by 
LIGO. 



In this section wc will provide rough estimates for determining the post-Newtonian orders, as a function of the 
(unequal) masses of the binary, necessary to construct templates that will be sufficient for precisely estimating the 
parameters of a detected source of gravitational waves. 

We begin by noting that the change in the number of cycles spent in a logarithmic interval of frequency can be 
computed within the post-Newtonian approximation and has the generic form |33j 



dN, 



eye 



r 



-5/3 



din/ 967r^5/3 



oo 
p,q=0 



(3.15) 



where / is the frequency of the emitted gravitational radiation, Cpq are coefficients that we assume to be 0(1) |47j . 
and ^ — G{m + Af )7r in units where c — I. The total number of cycles spent in a given frequency interval is found by 
integrating both sides of (3.151 so that 



5 / f 



p=0 



/ 



Cpq (q_5)/3 1 / „(g-5)/3 _ f{q-5)/3 

967r ^ g-5^ ' V''*^'=° ■' 

p,q=0 ^ 



p,q=0 



(3.16) 
(3.17) 



where we have used Kepler's Law = (,f to relate the (circular) relative velocity to the gravitational wave frequency 
and the prime on the second summation means that the g = 5 term is omitted. We take the upper limit of the 
frequency interval to be the gravitational wave frequency at the innermost stable circular orbit (ISCO). The ISCO 
demarcates the boundary between the quasi-circular inspiral and plunge phases of the binary's evolution. 



It is widely believed that (3.151 should be computed to 3.5PN order for the follow ing r eason [7]. Using Kepler's 
Law, the leading order term {p, q — 0) scales as so that the term in the sum of (3.151 that goes as (2.5PN) 
will give an 0(1) correction to Ncyc and, equivalently, the phase of h{t). However, an 0(1) correction implies that 
the template has "slipped" from the signal by an 0(1) number of cycles thereby reducing the chances of detecting the 
signal using matched filtering. Instead, one should compute Ncyc to 3PN, or better yet 3.5PN for higher precision. 
Furthermore, these higher orders vanish in the limit that v ^ implying that these corrections are potentially small. 



The p,q — term in (3.161 denotes the leading order (i.e., Newtonian) prediction for Ncyc and is often used for 



making rough estimates of the number of cycles spent in LIGO's bandwidth. The remaining terms are the post- 
Newtonian corrections from higher order gravitational, relativistic, spin, radiation reaction, and finite size effects. If 
we regard LIGO's lowest detectable frequency to be approximately 10 Hz |35] then the ratio 



i?(«)(10Hz)^ 



7vi^c'^(ioHz) 3e(«-^)/3 /L'^;"^" - io(«-^)/3 



7Vi^i-'^(10Hz) 



ln(/„eo/10) 



(3.18) 
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with q > 5 and Cpg = 0(1), gives an approximate measure of the corrections to the 0(1) number of cycles (i.e., the 
2.5PN correction) from higher PN terms. Note that the ratio is independent of p. 

From (3.181 we can estimate the approximate PN order, with a specified precision, that ( |3.16 1 should be expanded 



through. If S denotes this precision then the PN order is approximately given by the solution of 

Ri%{10Rz) = S (3.19) 

for q/2. A more accurate treatment should be given by more quantitative means [H [33] but since Ncyc is currently 
known only through 3.5PN [5', T] then we will continue with this somewhat crude estimation. 

Figures l(a)|l(c)" show contour plots of the solutions to (3.191 as a function of the binary's total mass and the 



mass ratio for S — 0.1,0.05, and 0.01, respectively. For a precision of w 10% Figure 1(a) indicates that 3.5PN 
expressions are adequate for equal mass binaries with a total mass as high as ~ 80Mq but that 4PN expressions may 
be more suitable for total masses between ~ 80Mq — 43OM0. Equal mass binaries with a total mass > 430Mq emit 
gravitational waves at the ISCO with frequencies less than lOHz and therefore do not enter LIGO's bandwidth. At 



a precision of 1% Figure 1(c) shows that for a binary with a mass ratio of 1/10 and a total mass of IOOOM0 that 
one may need PN expressions to order 6.5PN. For a precision of 10% Figure [T(a)| implies that the same system may 
require 4.5PN corrections. Such high PN orders for this particular binary should be expected since its waveform is in 
LIGO's observable bandwidth for a short frequency interval. 

The plots in Figure ([T|) should be regarded as providing crude estimates for determining the PN orders needed 
for LIGO to measure the number of cycles in its bandwidth with a given precision. Of course, these estimates can 
be justified only after the PN expansion of N^yc has been performed to higher orders in conjunction with a more 
comprehensive analysis. Nevertheless, we believe that Figure ([ij and the considerations of this section can serve as a 
motivation for calculating higher order corrections in the PN expansion for binaries with unequal masses. 

IV. FURTHER DEVELOPMENTS 
A. Gravitational perturbation theory with self-consistent backreaction 

The GP-SCB approach provides a general framework to build perturbation theories that are specific to the assump- 



tions and input being considered. We demonstrated in Section II B how the well-known post-Newtonian expansion 
and perturbation theory for EMRIs can be derived in the GP-SCB formalism. 

In Section [II C| we used the GP-SCB approach to introduce an adiabatic expansion with backreaction for arbitrary 
masses and velocities. Perhaps the most pertinent issue is determining the domain of validity of the equations of 



motion (2.57l-(2.62l. In particular, the degree of accuracy of those equations needs to be determined by the data 
analysis requirements set by gravitational wave interferometers. Can these system of equations adequately capture 
the plunge and possibly merger phases? How does the mass ratio and relative velocity aff'ect the answer to these 
questions, if at all? We hope to address these questions in some more detail in future papers [^ . 

A related issue is obtaining the solutions to such complicated backreaction equations. Since all of the degrees of 
freedom in the example of Section |II C[ are interacting non-trivially and non-linearly with other variables it would 
seem that there is a rich variety of processes that makes calculating the number of cycles that a binary stays in 
a detector's bandwidth, for example, difficult to perform and subsequently to compare with known results, both 
analytical and numerical. Furthermore, in many cases there are isometrics associated with the background geometry 
that can be exploited to assist in developing (semi-) analytical solutions. This is the case with the post- Newtonian 
method applied to the Minkowski spacetime, which is maximally symmetric, and with EMRI-perturbation theories 
off a Kerr background, which is axisymmetric. However, there does not seem to be any isometry associated with the 



background metric gofii^ since its evolution is determined in part by sources lacking any specific symmetry; see (2.571. 
Consequently, we will likely be forced to introduce some additional assumptions, but hopefully not as drastic as weak 
field, slow motion or small mass ratio, that could make it easier to find solutions beyond these simpler and more 
accessible regimes and provide new insights into processes hitherto too difficult to comprehend analytically. 

There have been several approaches that are based on an adiabatic approximation from a two-time scale separation 
for EMRI binaries [34, 35^^36, ,37, ,38, JE,. While our method is not restricted to extreme mass ratios, it is nevertheless 
more involved because the background geometry evolves with the backreaction from both the stress energy of the 
compact object (with finite size effects due to the extended nature of the mass) and the effective stress energy of the 
gravitational waves. 

The generic framework for building perturbative expansions, originally developed in |20j . may be useful for gener- 
ating other approximations in addition to the one discussed in Section [II C[ Depending on the dynamical regime of 
the two-body dynamics it may be possible to develop more perturbation theories suitable for other scenarios involving 
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gravitational binaries including the head-on collision of two black holes, the merger of comparable mass compact ob- 
jects, etc. One obvious perturbation theory that can be constructed would use the effective point particle description 
for both masses of the binary and derive the equations of motion in a similar manner as in Section [II C| Such a theory 
might lend itself useful to more practical schemes for solving the backreaction equations. We plan to investigate these 
and other issues associated with the self-consistent approach in future papers [22^ 23j . 



B. Post-Newtonian effective field tlieory for arbitrary mass ratios 



In Section III we provided crude estimates indicating the PN order that the number of gravitational wave cycles 
falling into LIGO's bandwidth should be calculated through for several given precisions. The PN order depends 
non-trivially on the binary's mass ratio and total mass. For a binary composed of IOMq and lOOM© compact objects. 
Figure 1(b) suggests that one may need to calculate N^yc to 4.5PN to achieve a ^ 5% precision relative to the 2.5PN 
correction. We also find that for a constant total mass the PN order increases as the mass ratio decreases. Meanwhile, 
for a constant mass ratio the PN order increases as the total mass increases. Both results are expected on intuitive 
grounds. 

Justifying these estimates requires computing to higher PN orders beyond the current 3.5PN expressions. In the PN- 
EFT approach, such calculations can be accomplished systematically and relatively efficiently compared to standard 
formalisms. Indeed, where higher PN orders are needed new physical interactions can manifest, including finite size 
effects from companion-induced tidal deformations (i.e., the non-effacement of internal structure), spin effects with 
radiation reaction, etc. We intend to investigate some of these features in future work. 
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